Anomalous Hall conductivity: local orbitals approach 



O 

(N 

< 



5-H 

CD 



C3 



I 

C 

O 

o 



> 

o 
in 

O 
O 



X 



Pavel Stfeda 

Institute of Physics, Academy of Sciences of the Czech Republic, 
Cukrovarnickd 10, 162 53 Praha, Czech Republic 
(Dated: April 12, 2010) 

A theory of the anomalous Hall conductivity based on the properties of single site orbitals is 
presented. Effect of the finite electron life time is modeled by energy fluctuations of atomic-like 
orbitals. Transition from the ideal Bloch system for which the conductivity is determined by the 
Berry phase curvatures to the case of strong disorder for which the conductivity becomes dependent 
on the relaxation time is analyzed. Presented tight-binding model gives by the unified way exper- 
imentally observed qualitative features of the anomalous conductivity in the so called good metal 
regime and that called as bad metal or hopping regime. 
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I. INTRODUCTION 

It has been known for more than a century that ferro- 
magnetic materials exhibit, in addition to the standard 
Hall effect when placed in a magnetic field, an extraor- 
dinary Hall effect which does not vanish at zero mag- 
netic field. The theory of this so-called anomalous Hall 
effect has a long and confusing history, with different ap- 
proaches giving in some cases conflicting results. While 
more recent calculations have somewhat unified the dif- 
ferent approaches and clarified the situation, it is still an 
active topic of research, see for example the recent review 
by Nagaosa et al^). 

It is generally accepted that anomalous Hall effect is 
induced by spin-orbit coupling. It was first suggested 
by Karplus and Luttinger 2 in 1954 to explain anoma- 
lous Hall effect observed on ferromagnetic crystals. Their 
analysis leads to the scattering independent off-diagonal 
components of the conductivity, which are attributed 
to the so-called "intrinsic" effect. Later, theories of 
this effect based on several specific models have been 
developed^. As shown recently, it is accompanied by a 
strong orbital Hall effed££. The conductivity is affected 
by scattering, which in the presence of spin-orbit cou- 
pling is basically of two types, the so called side-jump 7 
and skew scattering^—. They also lead to an anoma- 
lous Hall effect, called "extrinsic". It has also been ar- 
gued that in the anomalous Hall regime a periodic field of 
electric dipoles (electric polarizability) is induced by the 
applied current 2 - : 11 ' 12 . This property has recently been 
shown to be related to the so-called orbital polarization 
moment^ which is determined by Berry phase curvature 
in pure systems 1 ^. 

The best quantitative agreement with experimental ob- 
servations has been obtained by semi-classical transport 
theory 14 , leading to the Berry phase correction to the 
group velocity. For Fe crystals^ it yields an anomalous 
conductivity ~ 750 f2 _1 cm^ 1 while a value approaching 
1000 SI -1 cm -1 has been observed. However, up to now, 
generalization of this theory to systems with strong dis- 
order or subject to other types of fluctuations seemed to 
fail. It is the main aim of this paper to present theoretical 



treatment filling this gap. 

In contrast to standard transport theories the Hall con- 
ductivity is expressed in terms of local atomic-like or- 
bitals. It is explicitly derived that for perfect Bloch elec- 
tron systems this description coincides with that given 
by Berry phase curvatures indicating intrinsic, topolog- 
ical, origin of anomalous Hall effect. Presented view 
based on atomic-like orbitals allows to include effect of 
disorder by a local energy fluctuation of these orbitals 
which is an alternative description of scattering events. 
Fluctuations modify the Hall conductivity and transition 
from perfect to strongly disordered systems is analyzed. 
So called "intrinsic" and "extrinsic" Hall effect are just 
treated by an unified way. To test presented view based 
on local atomic-like orbitals the two band model within 
tight-binding approach is used. Obtained dependence 
of the anomalous Hall conductivity on the relaxation 
time shows the observed qualitative features 1 ^. Similar 
features has also been obtained by the different proce- 
dure with relaxation time being a fitting parameter 17 . 
In contrast to this work the presented model allows to 
relate relaxation time to fluctuations of atomic-like or- 
bitals, i.e. to its microscopic origin. By the unified way 
it gives experimentally observed scaling of anomalous 
Hall conductivity a xy with diagonal conductivity com- 
ponent <j xx in the region of the so called good metal for 
which <j xy ~ {<J XX )° and that of the bad metal for which 

The paper is organized as follows. In section II basic 
properties of electron systems in crystalline structures 
with spin-orbit coupling are summarized, while basic ex- 
pressions for Hall conductivity derived by using the quan- 
tum linear response theory are rederived in the section 
III. In section IV an alternative expression for Hall con- 
ductivity, including on-energy shell matrix elements only, 
is derived. Section V is devoted to perfect Bloch elec- 
tron systems where conductivity is expressed in terms of 
local orbital polarization moments which are further ex- 
pressed via the Berry phase curvature. In section VI a 
two-band model in tight-binding approach is presented 
to estimate the effect of a finite electron life time to the 
anomalous Hall conductivity at zero magnetic field. Pre- 
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sented theory of anomalous Hall effect is summarized and 
commented in the last section. 



II. SINGLE ELECTRON HAMILTONIAN AND 
STATISTICAL FLUCTUATIONS 

Within a mean field approach, electron properties are 
controlled by a single electron Hamiltonian which we con- 
sider in the following standard form 



H = 



V 



2mo 



\2 

V(f) + -4 a • W(r) x p 
in, 



[lb B eS ■ a . (1) 



Here, mo and \ib denote free electron mass and the Bohr 
magneton, respectively, p is the momentum operator, 
V{r) denotes a background potential and components 
of a are Pauli matrices. The third term on the right 
hand side represents spin-orbit coupling with A c being 
an effective Compton length. The last term on the right 
hand side describes Zeeman-like spin splitting due to the 
exchange-correlation energy represented by an effective 
field £> c ff which can generally be position-dependent. The 
corresponding velocity operator reads 



1 

ih 
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\2 

-| a x W(r) 



(2) 



Strictly speaking, Hamiltonian H quite well defines 
properties of electrons located within a finite volume of 
characteristic dimensions determined only by the elec- 
tron coherence length. Within each of such volumes the 
background potential as well B c g can be different. This 
way, a set of electron systems, the statistical ensemble, 
is defined. If time-dependent fluctuations can be treated 
within an adiabatic approach^, they can also be included 
in this ensemble. Measurable quantities are given by 
their statistically averaged values. It is useful to split 
Hamiltonian H, Eq. ([1}, into two parts 



H = H + SH 



H = (H) t 



(3) 



where only 5H depends on statistical fluctuations. For 
crystalline solids, to which the present treatment is de- 
voted, statistically averaged Hamiltonian H , obeys full 
crystal symmetry and the effective field B c ff can be as- 
sumed constant. It defines the so called virtual crystal 
with eigenstates \n, k) of the energy E n (k), characterized 
by the band index n and the wave vector k. Eigenfunc- 
tions are two-component Bloch spinors 



ikf 



H \n,k) = (H 



v 8tt 3 

\n,k) = E n (k) \n,k) 



(4) 



where u n Ar) are spinor functions periodic in r with pe- 
riod defined by the elementary lattice translations. Ve- 
locity matrix elements are diagonal in the wave vector k 



located within the Brillouin zone 

(n,k\v\n',k') = v n ' n '(k) S^p 



(5) 



and the expectation values read 



vjk) = (n,k\v\n,k) = - Vj;E n (k) . (6) 



Equilibrium properties are determined by the effective 
Hamiltonian, H e g, defined by the statistically averaged 
Green's function 



1 



H 



1 



z - H eS (z) 



(7) 



where z is the complex energy variable. It has the full 
crystal symmetry and it is diagonal in the representa- 
tion given by eigenstates of the averaged Hamiltonian 
Hq, Eq. (j4j. Effective Hamiltonian is non-Hermitian and 
energy dependent but it is analytic in both complex half- 
planes, H e g(z*) — H^~ s (z). Its standard form reads 

H eS (z) = Ho + E(«) , E(z) = A(z) - iT(z) , (8) 

where is the energy dependent self-energy. Inverse 
value of its imaginary part represents a finite electron 
life-time. 

To include an external magnetic field B the Hamilto- 
nian defined by Eq. ([1} has to be modified. Both the 
momentum operator entering the Hamiltonian, and the 
velocity operator, Eq. ([2]), have to be replaced by their 
counterparts, which include a vector potential A. 



V 



P 



i A . 



(9) 



Here, e denotes the electron charge absolute value and the 
magnetic field is given as B — curl A. Also the value of 
the parameter B c g defining Zeeman-like splitting is mod- 
ified by B. The external magnetic field generally removes 
translation symmetry. Exceptions are the so called ratio- 
nal magnetic fields for which the problem becomes invari- 
ant under translations with different elementary transla- 
tions than those given by the periodic potential. 



III. LINEAR RESPONSE THEORY 

In this section the standard linear response theory will 
be described and well known general formulae derived to 
summarize basic theoretical assumptions. In accord with 
the original work by Kubo^ the external electric field is 
supposed to be turned on at the time t — > — oo and reach 
the final value £ at t = 0. Exponential time evolution is 
considered, £(t) = £ e 6 *^, with e > being an infinites- 
imally small quantity. It gives rise to the Hamiltonian 
perturbation er ■ £ (t) and corresponding density matrix 
pit) has to satisfy the equation of motion which reads 



dpjt) 
dt 



1 

ih 



pit), H + ef-E(t) 







(10) 
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It can be written in the following way 



p(t) = po(H) + e 



-iHt/h 



Pl(0 



jHt/h 



(11) 



where pi (t) represents the deviation from the density ma- 
trix in the absence of electric field given by the Fermi- 
Dirac distribution function 



Po(H) 



H-p. \ — 1 

1 + e 'b t 



(12) 



where /i and fceT are the chemical potential and thermal 
energy, respectively. Up to linear terms in the electric 
field, pi(t) reaches the following value at t = 



+oo 



Pi(0) = ^ / < 



-et/% e —iHt/% 



Po(H) , ef ■ £ 



and the resulting current density reads 
j = — e lim (Tr[pi(0)wl 

€->0+ 



(13) 
(14) 



Generally, the above potential perturbation e£f of the 
Hamiltonian is not the only perturbation caused by elec- 
tric field. The potential gradient e£ enters the spin-orbit 
term of the Hamiltonian, Eq. ([T} , as well as the velocity 
operator, Eq. ([2J, giving rise to additive terms linear in 
electric field. Resulting contributions to the current den- 
sity are of a higher order than A 2 in the Compton length. 
For the considered Hamiltonian H, Eq. ((TJ), which in- 
cludes spin-orbit coupling only approximately up to the 
order of A 2 , these contributions thus have to be ignored. 

Introducing the (5-function operator 



S(rj - H) 



lim 



2m 



(15) 



1 



rj — H ± it 



the time integration, Eq. (JT5]), can easily be performed. 
For the limiting case of the zero temperature in function 
Po(H) the components of the conductivity tensor obey 
the following form 



coherent system is defined by the physically acceptable 
limit r — > + . That is why it is sufficient to consider just 
one infinitesimally small parameter e — > + implicitly en- 
tering Eq. (TT5|) . Misunderstanding of this limit has been 
reason for doubts concerning validity of resulting expres- 
sions, especially of the Hall conductivity. 

The non-zero temperature smearing effect of the dis- 
tribution function yields 



drj 



a tj (r])dr) . 



(17) 



This relation allows to limit our attention to the analysis 
of the energy-dependent conductivity given by Eq. (|16p . 
It does not mean that cry (/_t) has to be temperature inde- 
pendent since potential fluctuations as well as some of the 
Hamiltonian parameters can be implicitly temperature- 
dependent quantities, as is e. g. B c g representing the 
exchange-interaction effect. 

Using definition of the velocity operator 



1 

v = T 
in 

and the identity 



[r,H] 



i 

ih 



dG ± 
drj 



(18) 



(19) 



the Eq. (fTB)) can be rewritten^ in a form including ve- 
locity operators only 



<jij (n) = the x 
Tr I 5(r] - H) 



dG + 



dG~ 

Vi— — Vi 



drj 



drj .(20) 



It coincides with that obtained by Bastin et al2i by the 
use of a different procedure. For diagonal components of 
the conductivity tensor the integration by parts gives the 
well known Kubo-Greenwood formula 2 ^ 



au(n) = ire 2 h(Tr{v l 6(p- H)vi8{pL- H)}) B 
Making use of the following identit y 20 ' 23 



(21) 



e 2 / (Tr {,5(77 - H) [v.G+r, + r^Vi] }) av dr, , (16) 



where i,j = x, y, z. The proper way to treat the limits 
in Eqs. (fT4|) (due to the electric field time evolution) and 
([T5jl (regularization of 5(H — r])) would be to introduce 
two different infinitesimal parameters and treat both of 
them independently, after all other steps are taken. How- 
ever, our procedure of statistical averaging yields a non- 
zero imaginary part of the self-energy, T, entering both 
averaged operators, and the limiting case of the fully 



iftTr \ 5(r) - H) 
ih d 



dG+ 

Vi— Vj 

dr] 



dG~ 
drj 



= y ^ Tr {S(r, - H) [viG+^Vj - VjG-fj})^ } + 

+ l^{S(v-H)[r i v j -r j v i ]} ,(22) 

the off-diagonal conductivity components can be split 
into two parts^ 



(ii) 



(23) 
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where 



vy vMj — e ~2 

x (Tr {<S(,i - H) [v x G+(fi)v y - v y G-(»)v x ] }) av , (24) 



and 



^(M) = y (Tr {5(77 - if) [xvy - yv x ]}) av = 

= -e 2 (Tr{S( V -H)yv x }) av . (25) 

These conductivity formulae are quite general, they also 
include both the effect of the magnetic field and the 
spin-orbit coupling 25 . Similar treatment, but formally 
more complicated^, is applicable to other transport co- 
efficients, like the thermopower and the heat conductiv- 
ity. 

The above expressions for the Hall conductivity a xy (p) 
have been particularly useful for understanding of the 
quantum Hall effect. While the contribution cr X y(fj,) van- 
ishes whenever /i is located within an energy gap, the con- 
tribution <j X y* (p) can remain finite, giving rise to quan- 
tized values of the Hall conductivity 24 . However, these 
formulae are not very convenient for fj, located within an 
energy band since, at least for some specific model o 17 ' 26 , 
a non-negligible part of <T X y (m) i s nearly canceled by part 
of the a X y (p) . Thus it should be useful to derive an al- 
ternative form of the expression for the Hall conductivity, 
which should also allow a clear physical interpretation. 



IV. HALL CONDUCTIVITY IN TERMS OF 
ON-ENERGY-SHELL MATRIX ELEMENTS 

Conductivity can directly be measured on samples 
with so called Corbino disc geometry. In the limiting 
case of the large disc radii such samples can be approx- 
imated by strips with a rectangular cross-section L y L z . 
Using a proper choice of the vector potential, periodic 
boundary conditions along x direction can be considered 
on the length L x . This geometry allows to apply an elec- 
tric field along y direction, and establish <J yy , and also the 
Hall current along y direction and consequently a xy . Un- 
der these boundary conditions the eigenvalue problem for 
one particular Hamiltonian of the form given by Eq. ([T]) 
can be solved, at least in principle. Obtained eigenstates 
| a) of energy E a represent one of the systems belonging 
to the considered statistical ensemble. A corresponding 
contribution to the Hall conductivity can be analyzed by 
using this a-representation. 

The term given by Eq. (|24p reads 



(a,E) 



ih \ -> 
~2 ^ 



5(E - E a 



where u">™ denotes velocity matrix elements. Inserting 
commutation relation Eq. (|I8p for the operator v y , we 
get 

<4?(«> E) = S(E - E a ) (E a - E a .) x 



a, a' 

x y 

E - E n i + i 



„,a,a a ,a 
y u x 

E-E„,-i 



(27) 



Note that because of the considered strip geometry, ma- 
trix elements of the y coordinate are finite. Since terms 
for which E a — E a i vanish, we get 



a ( i y \a,E)= e -Y,KE-E a ) 



v y +y v x 



i 



J E a ,,E a 



(28) 



In a-representation the remaining conductivity contribu- 
tion given by Eq. (|25|) can be written as 



(29) 



ag I \a,E) = ~ J2$(E-E a ) 



a, a a .a , a, a a ,a 

v x ' y +y v x > 



(26) 



Sum of cr^y an d v X y yields the Hall conductivity ex- 
pressed in terms of on-energy-shell matrix elements only 

a xy (E) = -e 2 /j2 5 ( E - E ^ < a 'y a '' a 5 B a ,,E a \ , 

\a,a> / av 

(30) 

where averaging is taken over the a-representations of all 
elements of the considered ensemble. In the case that the 
electron system is fully coherent within the considered 
sample volume L x L y L z the averaging procedure should 
be avoided. 

The above expression for the Hall conductivity, 
Eq. ([30| . has clear interpretation for quantum Hall effect. 
Let us assume that the Fermi energy is located within re- 
gion of localized states. From their definition only diag- 
onal matrix elements of the y coordinate can differ from 
zero, y a ,Q = y a,a 5 a ' t0l . Since their velocity expectation 
values vanish, contribution of localized states to the Hall 
conductivity vanishes as has to be. Non-zero contribu- 
tion can only be given by chiral edge states giving rise 
to quantum Hall effect non-affected by the presence of 
localized states. 



V. PERFECT BLOCH ELECTRON SYSTEMS 

It is of particular interest to apply expression for the 
Hall conductivity, Eq. (PU)) , to systems having translation 
symmetry, including those modified by rational magnetic 
fields. The basic property of these systems is that the 
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velocity matrix elements are diagonal in wave vector k, 
Eq. ([5]). In general, the probability to find on the Fermi 
surface states having the same k but belonging to differ- 
ent bands is statistically negligible, except of the trivial 
case of band degeneracy. Since Eq. (|30|) includes only 
on-energy-shell matrix elements, inter-band matrix ele- 
ments do not affect conductivity of ideal Bloch electron 
systems. 

To proceed further, let us first discuss the contribution 
to the Hall conductivity, Eq. ([50]). of states having zero 
velocity expectation value along y direction. It is given 
by 

eV ws (v^ n (k)y n (k) + v ^ n (-k)y n (-k)) , (31) 

where V ws denotes volume of the Wiegner-Seitz cell. It 
is a contribution of the local orbital momentum to the 
conductivity. As has been already discusse d 13 ' 27 it is re- 
sponsible for charge polarization in transport regime and 
that is why it can be called orbital polarization momen- 
tum. States with non-zero velocity along y direction con- 
tribute by the product - 
can be written as 



and the Hall conductivity 



(E) 



(32) 



where the orbital polarization momentum P n {E) reads 
—7T~~~ = - A / S (E n (k) - E) f n (k) x v n (k ) d 3 k . 

BZ 

(33) 

Here r n (k) denotes the expectation value of the radius- 
vector of the orbital n within the Wiegner-Seitz cell and 
integration is limited to the Brillouin zone volume. 

It can be easily shown that P n (E) can be expressed in 
terms of the Berry phase curvature. Since 

6{E - E n (k) ) = , (34) 

dE n (k) 

where fo(E) denotes the zero-temperature Fermi-Dirac 
distribution, we get 

P n(E) =-^~ I [Vg/o (£„(£))! xr n (k)d 3 k, (35) 



hAir 2 



BZ 



Integration by parts gives 
P n (E) = J f (E n (k)) (Vg x f n {k)) d 3 k. (36) 



BZ 



Using the expression derived in the Appendix for expec- 
tation values of the radius- vector 



.(*) 



Im 



f/ 3 7 



(37) 



we finally get 



BZ 



where O n (fc) is just the Berry phase curvature^ 



Jln(fc) = -Im / [Vgti+^J x (Vg^f)) d*r . (39) 



To conclude, the effect Berry phase curvatures is an al- 
ternative description to the presented effect of the orbital 
polarization moment, Eq. (|33|) . 



VI. 



TIGHT-BINDING APPROACH AND 
SINGLE SITE FLUCTUATIONS 



The aim of this section is to present a simple model 
system allowing us to understand the main features of 
the anomalous Hall conductivity at zero external mag- 
netic field. For the sake of simplicity the consideration 
will be limited to isotropic systems, like those of the cu- 
bic symmetry. A single band model Hamiltonian will be 
considered of the following form 



H = J2V) (E a +Si){l\ 



(40) 



l.rn 



where \l) and |m) are Wannier functions representing 
atomic-like orbitals associated with lattice sites Ri and 
R m , respectively. To model fluctuations, a variation 5i of 
atomic- like orbital energies, E ai will be considered, while 
hopping integrals ti, n are supposed to be fluctuation in- 
dependent quantities. Instead of considering a specific 
form of the band dispersion for averaged Hamiltonian, 
Ha = (-ff)av = H{5 — > 0), we shall assume it gives an 
elliptical density of states. In accord with Hubbard 2 ^, its 
energy dependence normalized per Wiegner-Seitz volume 
Vws can be written as follows 



9(E) = 
9(E) 



2 

7TW 2 



\E-E a \>w, 
y/w 2 -(E- E a ) 2 , \E~E a \<w, 



(41) 



where 2w denotes the band width. The corresponding 
mean Fermi velocity reads 



Up 



w 
2h 



{\ ) 2/3 



(42) 



1 /3 

where a = Vws just equals to the lattice constant for the 
simple cubic lattice. Note that this model system was 
already successfully applied in a description of electron 
properties of substitutional alloys 29 . Under these sim- 
plifying assumptions the single site orbital polarization 
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moment denned by Eq. (|33[) becomes an energy depen- 
dent quantity, given as 



lattice site, we get 



Pi{E)=n a ~a 2 - 



(f) 1 



r a (Si) 



. (43) 

where r a (Si ) denotes the average radius of the considered 
atomic- like orbital of the energy E a +Si and n a is the unit 
vector parallel to its orbital momentum. The only free 
model parameter in Pi(E), and consequently in anoma- 
lous Hall conductivity, Eq. (j3"2")l , is the orbital radius. For 
a ~ 3A and r a (0)/a ~ 0.3, the corresponding anomalous 
Hall conductivity can reach values of several hundred of 
f2 _1 cm _1 , similar to those observed experimentally. 

In real structures several overlapping energy bands 
contribute to the conductivity. Let us for simplicity con- 
sider two bands of the same width originated in atomic- 
like orbitals of energies Eq ± AE a . To model ferromag- 
netic state we assume that B c g = (0, 0,-B G ff), and that 
electron states belonging to different bands have opposite 
spin orientation along z direction, and opposite orienta- 
tion of their local orbital moments. With rising energy of 
atomic orbitals their radius increases. Up to the lowest 
order in energy we get 



r a (E a ) 



= 2ll 



E a — E 



w 



(44) 



where r denotes the radius for the orbital state of the 
energy Eq and the parameter k represents how the radius 
changes with the orbital energy. 

The origin of fluctuations modelled by local energy 
shifts S can vary considerably. In addition to the con- 
figurational disorder like impurities, alloy composition 
and other types of the disorder there are thermal fluctu- 
ations and that given by local values of B e g. Assuming 
that electron hopping between lattice sites is indepen- 
dent of their band origin, fluctuations of orbital energies 
Eq ± AE a can be treated as independent. Ensemble av- 
eraging can thus be performed separately for each of the 
band contributions. To get a qualitative estimation of 
the effect of fluctuations the Lorentzian distribution of 
the parameter d will be considered 



1 



7T s 2 + r 2 ' 



(45) 



where the energy parameter Y defines strength of the 
fluctuations. For the sake of simplicity so called virtual 
crystal approximation 2 ^ is used to establish conductiv- 
ity. In such a case the parameter T represents imaginary 
part of the self-energy £ entering effective Hamiltonian, 
Eq. ||5J|, while the real part of the self-energy equals to 
zero. It also defines the electron life time, r = H/2Y. 
Substituting Eq. (|4"3"f into Eq. ([3^)1 for the anomalous 
Hall conductivity, and assuming there is one electron per 



(T xy (n = E ) = — - — x 



ah V2 



2/3 



+1 



x / (1 - e 2 )(l + <) b(£ + e,7) - P« - e, 7)] d£ , (46) 



where the half band width w has been used as the energy 
unit, 7 = Y/w, e = AE a /w. The resulting conductivity 
dependence on the fluctuation strength represented by 
the parameter 7 = Y/w is shown on Fig. [T] for the fol- 
lowing set of parameters: a — 3A , r a (0)/a — 0.3, K = 1, 
and e = AE a /w = 0.6. 




10' Y=r/w 10° 



FIG. 1: Anomalous Hall conductivity for two band model as 
a function of the fluctuation strength F = h/2r. Inset shows 
corresponding density of states. 

The anomalous Hall conductivity dependence on the 
fluctuation strength T shows the same qualitative fea- 
tures as that found for multi-d-orbital tight-binding 
model developed by Kontani et alii. In the case of a 
weak disorder, T/AE a <C 1, conductivity is nearly con- 
stant while for strong disorder, T / AE a > 1, it decreases 
with a power of the electron life time r = h/Y, which 
for the present example even slightly exceeds quadratic 
dependence. In comparison with the procedure used by 
Kontani et al whose evaluation is based on the conductiv- 
ity formula including only velocity operators, presented 
treatment based on local orbital polarization moments is 
much simpler and has a clear interpretation. Ferromag- 
netic state, necessary for appearance of anomalous Hall 
conductivity, can be characterized by average spin along 
magnetization axis which is supposed to be parallel with 
£ direction. In the case of the considered two-band model 
it is assumed that each of them is fully spin-polarized but 
in opposite directions, s z = ±1/2. If there are no fluctu- 
ations the average spin per site (s z ) reaches a maximum 
value. Fluctuations of orbital energies can only lead to 
suppression of this value. The dependence of (s z ) on 
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the parameter T representing fluctuation strength is pre- 
sented in Fig. [5] for the same set of parameters as that 
used for dependence presented in Fig. Q] It shows the 
same qualitative features as the dependence of the Hall 
conductivity on T. In the region where (s z ) is smaller 
than its maximum value the Hall conductivity depends 
on the fluctuation strength. It is the region in which the 
disorder is strong enough to weaken ferromagnetic state. 
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FIG. 3: Scaling of the Hall conductivity a xy with longitudinal 
conductivity a xx . Full line is given by the same model param- 
eters as that used for the dependence presented in Fig. [1] while 
for dashed line the parameters n and e = AEa/w have been 
changed (« = 0.5 and e = 0.4). Dashed-dotted line shows the 
slope 1.6. 



FIG. 2: Average spin (s z ) per site for two band model as a 
function of the fluctuation strength V — h/2r. Model param- 
eters are the same as that used for the dependence presented 
in Fig. [[] 

The simplest approach for diagonal conductivity com- 
ponent <r xx , Eq. (1211) . is to neglect vertex corrections, i.e. 
to use decoupling (GvG) av — > (G) av v(G)i lv . In this ap- 
proximation relaxation time coincides with electron life 
time r = H/2T. Using virtual crystal approximation for 
the considered two-band model and probability distribu- 
tion of fluctuations, Eq. (|21[) . evaluation of a xx becomes 
trivial. In Fig. [3] obtained scaling of the Hall conductiv- 
ity with a xx is shown. It reveals typical features observed 
experimentally^^. Especially the case of the moderate 
disorder (good metal regime) for which \a xy \ ~ (<J XX )° 
and that of the strong disorder (bad metal / hopping 
regime) for which \a xy \ — (a xx ) 1,6 . 



VII. CONCLUDING REMARKS 

Transport theories of conductivity arc traditionally for- 
mulated to give expressions containing velocity matrix el- 
ements only. In the presented treatment the Hall conduc- 
tivity has been expressed in terms containing also posi- 
tion operator. It has been found that only on-energy-shell 
matrix elements of both operators are relevant, Eq. (1301) . 
For Bloch electron systems the Hall conductivity is given 
by a part of the orbital magnetization of Fermi electrons, 
called orbital polarizability momentum, Eqs. (|32p and 
(I33p . which is a quantity determined by atomic-like or- 
bitals. For a perfect Bloch system it is equivalent to the 
expression given by the Berry phase curvatures, Eq. (13"8"1) . 

To test applicability of the derived alternative expres- 
sion for the Hall conductivity, Eq. (|30p. the two-band 
model based on the tight-binding approach has been 



used. Disorder has been modelled by energy fluctuations 
of single site orbitals. It represents intra-band scatter- 
ing which can be identified with the so called side-jump 
scattering. However, this simple model excludes effect of 
the skew scattering since spin of electrons is supposed to 
be fixed. Despite of its simplicity it correctly describes 
scaling of the anomalous Hall conductivity a xy with diag- 
onal conductivity component a xx in the region covering 
the so called good metal regime, \cr xy \ ~ (cr xx )°, and 
the bad metal (hopping) regime, \a xy \ ~ {o xx ) l ' e ■ The 
bad metal regime has been identified with the regime in 
which disorder becomes strong enough to weaken ferro- 
magnetic state. It remains an open question if the pre- 
sented form of the general expression for the Hall conduc- 
tivity in terms of the local orbitals could also be effective 
in description of the skew scattering effect. 

To estimate anomalous Hall conductivity for a real ma- 
terial requires knowledge of local orbitals represented by 
Wannier functions and also the specific form of fluctu- 
ations relevant for the studied system. In particular, 
finite temperature spin fluctuations are expected to be 
essential. It is a challenge to work out such a procedure 
based on the first principle calculations. Newly developed 
numerical techniques allowing to establish Wannier func- 
tions giving the best tight-binding model parameter a 30 i 31 
seem to be a proper way to establish the fluctuation effect 
upon the anomalous Hall effect in real materials. 
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Appendix A 

Periodic part of Bloch functions can be expressed via 
Wannier functions <j> n (f— Ri) as follows 



1=1 



where by definition Wannier functions are orthonormal 
with respect to their mass-center position vector Ri . The 
expectation value of the position vector r reads 



T n (k) = (n,k\rn,k) = 
— [ \u + ~(r)ru ?(r) + u Ar) ru + Jr) 

2 / L ™, fe n ' k ™' fe »».* 



dr 6 = 



N 



1=1 
N 



i=i 



U r (r) =- e 

* r-Ri 
2VN 



ERi 
9~/V 



27V 



ik(R l -R l ,)i + (i f_ 



ct>+(r-R l ,)(t> n {r-R l )d iS r+ 



N 

E 



1 N 



where the last constant term represents the center of 
mass of the considered electron system. For states with 
— k the above relation gives 



r i N 

f n (-k) = Im u+ u n> _ g (»0 d 3 r + - E ^ • 

In the summation over Fermi surface states of the prod- 
uct f n (k) v n (k) entering Eq. (|33p the mass-center posi- 
tion of the electron system is cancelled out since v n (k) = 
—v n (—k). The last term on the right-hand side of the 
above equation can thus be excluded from the consider- 
ation since it docs not affect the final result. 
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